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ABSTRACT 

Host star metallicities have been used to infer observational constraints on planet formation through¬ 
out the history of the exoplanet field. The giant planet metallicity correlation has now been widely ac¬ 
cepted, but questions remain as to whether the metallicity correlation extends to the small terrestrial¬ 
sized planets. Here, we report metallicities for a sample of 518 stars in the Kepler field that have no 
detected transiting planets and compare their metallicity distribution to a sample of stars that hosts 
small planets ( R p < 1.7 R®). Importantly, both samples have been analyzed in a homogeneous manner 
using the same set of tools (Stellar Parameters Classification tool; SPC). We find the average metallic¬ 
ity of the sample of stars without detected transiting planets to be [m/H] SNTP dwarf = — 0.02±0.02 dex 
and the sample of stars hosting small planets to be [m/H] STp = —0.02 ± 0.02 dex. The average 
metallicities of the two samples are indistinguishable within the uncertainties, and the two-sample 
Kolmogorov-Smirnov test yields a p-value of 0.68 (0.41 a), indicating a failure to reject the null hy¬ 
pothesis that the two samples are drawn from the same parent population. We conclude that the 
homogeneous analysis of the data presented here support the hypothesis that stars hosting small 
planets have a metallicity similar to stars with no known transiting planets in the same area of the 
sky. 

Subject headings: planetary systems — techniques: spectroscopic, spectroscopic surveys 


1. INTRODUCTION 

The discovery of thousands of transiting exoplanets 
has enabled studying the properties of statistically sig¬ 
nificant ensembles of planets and their host stars. Such 
studies, accompanied by planet formation theories, can 
be used to shed light on how planets form. In early 
stages of the exoplanet held when only a handful of plan¬ 
ets were known to exist, a correlation between host star 
metallicity and the prese nce of hot-Jupiter type planets 
was suggested (!Gonzalezlll997l) . This tendency was later 
confirmed by a series of papers showing that metal-rich 
stars are much more likely to harbor gas- giant planets 
(jSantos et al.ll200il : I Fischer fc Valentill2005t h This result 
is usually interpreted as natural support for the core ac¬ 
cretion model, where planets growing in a metal-rich en¬ 
vironment will have a higher likelihood of reaching the 
critical core mass allowing run-away gas accreti on before 
t he gas in the system has dissipated (|Ida fc Linl 120051: 
iMordasini et all[200 9. 2012). Recent discoveries of thou¬ 
sands of small exoplanets have allowed these types of 
studies to extend into the terrestrial planets regime. 

The hot-Jupiter type exoplanets with periods shorter 
than about 10 days dominated the initial discover¬ 
ies, but we now know these types of planets to 
be rare (found around ~ 1% of Sun-like stars; 
IWright et all 12012 1. On the other hand, the Ke¬ 
pler Mission has demonstrated that small planets are 
astonishingly common and ubiquitous in our galaxy 
([Howard et al.ll2012t iFressin et al.1120131 iPetigura et all 
120131 : iDressing fc Charbonneaul 1 201. 'll) thus supporting 
and extending the earlier hints from radial velocity sur¬ 
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veys t hat low-mass planets are common dMavor fc Udrvl 
120081 : iHoward et al.l l200l : iBonfils et al.1120131) Initially, 
the metallicity correlation was assumed to pertain to ex¬ 
oplanets in general. However, a number of later studies 
have shown that the giant-planet metallicity correlation 
does not seem to e xtend into the small p l anet regime 
dSousa et all 120111: IBuchhave et all l2012t iNeves et al.l 
120131 : lEverett et al M2013T) . Recent results measuring the 
metallicities of a large number of stars hosting small plan¬ 
ets suggest that the planets can be divided into three 
regimes, namely terrestrial planets, gas-dwarf planets 
small planets with low er mean densities) and gas-giants 
Buchhave et al. I2014D . where the planets in the terres¬ 
trial planet regime have a metal licity consistent with so - 
lar. Then, on the other hand, I Wang fc Fisclierl (|2015H 
used photometricall y derived metallicities from the Ke¬ 
pler Input Catalog dBrown et al.l 120111) com bined with 
spectroscopic metallicities from IBuchhave et al.l (12014H 
to suggest that the occurrence rate of small planets is 
higher around metal rich stars, thus claiming that the 
planet-metallicity correlation does extend all the way 
down to the small planets. 

The giant planet-metallicity correlation was deter¬ 
mined by examining the metallicities of a sample of host 
stars harboring hot-Jupiters and comparing it to a “con¬ 
trol sample” of stars that did not host such planets. This 
made sense, since the Doppler surveys could readily de¬ 
tect these massive close-in companions except for the rare 
cases where the orbital plane of the system was nearly 
face-on. A control sample of stars without hot-Jupiter 
companions was thus easy to construct for a volume lim¬ 
ited survey. As the radial velocity surveys increased their 
sensitivity and were able to detect sub-Neptune sized 
planets, it was still possible to construct a sample of stars 
that did not host planets of the detected type and esti¬ 
mate how the planet occurrence rate was affected by host 
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star metallicity (e.g. iSousa et al.ll 201 lD . 

However, it is much more difficult to construct a con¬ 
trol sample for the small transiting exoplanets discovered 
by Kepler. We know that small planets are common 
around solar-type stars, and while the transit method 
has the required sensitivity to detect such small planets, 
the orbital plane must be oriented such that the planets 
transit for us to detect them. As a result, we cannot con¬ 
struct a sample of stars free of planets from using transit 
surveys. The current level of precision of the leading 
radial velocity spectrographs is not yet able to detect 
the small masses of these planets, except in rare circum¬ 
stances. In fact, it turns out to be much more difficult 
to construct a sample of stars where we can definitively 
rule out the presence of small planets than to discover 
the planets themselves. 

In this paper, we examine a sample of stars in the Ke¬ 
pler field with no detected transiting planets (Stars with 
No detected Transiting Planets: SNTP) and compare 
their metallicities to stars hosting planets with sizes be¬ 
low R p < 1.7 f ?0 (Stars with detected Transiting Planets: 
STP). 

2. SAMPLES AND OBSERVATIONS 

Our sample of SNTP consists of 51 8 stars pre viously 
analyzed with asteroseismology ([Chaplin et all 12014 ). 
The stars were selected for observation by the Kepler 
Follow-up Program team (KFOP) as a sample of “gold 
standard” stars that would be well-characterized by both 
asteroseismology and spectroscopy. The stars were ob¬ 
served with the Tillinghast Reflector Echelle Spectro¬ 
graph (TRES) on the 1.5 m Tillinghast Reflector at the 
Fred Lawrence Whipple Observatory on Mt. Hopkins, 
Arizona using the medium resolution fiber (2 , /3 projected 
diameter) with a resolving power of R ~ 44 000, giving a 
wavelength coverage of ~ 3800 — 9100 A. 

Most of the stars were observed only once, but in some 
cases, several observations were available. In total, we 
analyzed 834 spectra of 518 stars. Since the magnitude 
of the host stars ranged from 3rd to 13th magnitude, the 
exposure time ranged from 45 seconds to 1 hour with a 
mean exposure time of approximately 11 minutes. The 
resulting average signal-to-noise ratio per resolution ele¬ 
ment (SNRe) in the Mg b region was 49 corresponding 
approximately to a signal-to-noise ratio per pixel (SNR) 
of 30. The spectra used in this paper are available on 
the Kepler Community Follow-up Observing Program 
(CFOP) website (http://cfop.ipac.caltech.edu). 

3. STELLAR PARAMETERS 

The sample of SNTP stars was analyzed with the Stel¬ 
lar Parameters Classificati on tool (SPC) i n the same 
manner as document ed in iBuchhave et, all (|2012j ) and 
iBuchhave et all (120141) . Here, we provide a brief overview 
of SPC. 

SPC cross-correlates the observed spectrum against 
a grid of synthetic template model spectra with vary¬ 
ing values of stellar parameters for effective tempera¬ 
ture, surface gravity, metallicity and rotational veloc¬ 
ity. The normalized cross-correlation function (CCF) 
peak height is determined for each template, indicating 
how well the observed spectrum matches the various syn¬ 
thetic model spectra. The library consi sts of an exten¬ 
sive grid of stellar models (iKuruc/JI 19921) and covers the 
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Figure 1. The metallicities of the stars in the SNTP sample de¬ 
termined using SPC with all parameters free versus SPC using 
the surface gravity from asteroseismology as a prior. The mean 
difference is -0.01 dex and the RMS is 0.07 dex, suggesting that 
metallicities can be determined reliably without a prior on surface 
gravity. 

wavelength region between 5050 — 5360 A. It spans an 
interval for effective temperature of 3500 K < T e ff < 
9750 K, surface gravity of 0.0 < log(g) < 5.0, metal¬ 
licity of —2.5 < [m/H] < +0.5 and rotational velocity of 
0 km s^ 1 < Kot < 200 km s -1 and has a spacing of 250 K 
in effective temperature, 0.5 in log(g), 0.5 dex in [m/H] 
and progressive spacing in rotational velocity yielding a 
grid consisting of 51,359 synthetic spectra. 

Rather than selecting the best matched template, 
which would restrict the stellar parameters to the rather 
coarse grid spacing, SPC assumes that the normalized 
CCF peak height varies smoothly between grid points. 
The CCF peaks are fitted with a three dimensional third 
order polynomial as a function of effective temperature, 
surface gravity and metallicity, weighting the peaks pro¬ 
portional to their height, giving less weight to the CCF 
peaks that are a poor match to the spectrum. The peak 
of the fitted surface is determined, yielding the final stel¬ 
lar parameters. 

SPC utilizes the entire spectrum in the wavelength 
region of the library, enabling the technique to derive 
stellar parameters from spectra with relatively low SNR 
(down to a SNRe of roughly 30 corresponding to a SNR 
of 18 per pixel). Since SPC takes advantage of all the 
absorption lines in the wavelength region, we denote the 
metallicities in this work by [m/H], representing a mix of 
metals assumed to be the same as the relative pattern of 
the abundances in the Sun, not to be confused with the 
abundances for individual elements (e.g. [Fe/H]). 

SPC determines all the stellar parameters simultane¬ 
ously, however, if some parameters are available from 
a more reliable source, we can use this information to 
set a prior on these parameters in SPC. The surface 
gravity is notoriously difficult to determine spectroscop¬ 
ically (e.g. [Torres et al.ll2012l) . so we used t he value de¬ 
termined from asteroseismology in [Chaplin et all (|2014f) 
as a prior. Since the sample of STP stars does not have 
surface gravities determined by asteroseismology, we in¬ 
vestigated whether enforcing a prior on the surface grav¬ 
ity would significantly affect the metallicities by running 
SPC with all parameters allowed to float and using the 
asteroseismic surface gravities as a prior. The resulting 
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Table 1 

Metallicities of star with no detected transiting 
planets 


KIC 

[m/H] 

1430163 

-0.05 ±0.08 

1435467 

0.01 ±0.08 

1725815 

-0.07 ±0.08 

2010607 

-0.01 ±0.08 

2306756 

0.42 ± 0.08 


Note. — Table [I] is published in its entirety in the 
electronic edition of the Astrophysical Journal. 
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Figure 2. The distribution of metallicities of the stars with no 
detected transiting planets (blue, 518 stars ) and the sample of star s 
with small planets (red, 157 stars) from iBuchhave et al 3 {mi). 
The vertical lines show the average metallicity of the two samples: 
[m/H] gNTp = —0.02±0.01 dex and [m/H] gTp = —0.02±0.02 dex, 
where the uncertainty is the standard error of the samples. 

metallicities are shown in Figure |T| The mean metal¬ 
licity difference between the two sets of metallicities is 
-0.01 dex and the RMS is 0.07 dex, indicating that us¬ 
ing a prior on the surface gravity has little effect on the 
metallicities. The metallicities of the SNTP stars are 
listed in Table [Q 

In order to determine the mass and radius of a star, 
asteroseismology relies on effective temperatures and 
meta llicit ies f rom e.g . spectroscopy or photometry. In 
IChaplin et al.l (12014I) . photometrically derived effective 
temperatures were available for the entire sample, but 
since no reliable metallicities were available, the metallic- 
ities were set to —0.2 ± 0.3 dex. Spectroscopic estimates 
of effective temperature and metallicity were available 
for a subset of 87 stars. In this paper, we focus on the 
metallicities of the host stars, but a second iteration of 
the asteroseismic results using stellar parameters from 
SPC with a asteroseismic prior on the surface gravity is 
the subject of a paper in preparation where all the stel¬ 
lar parameters will be published. However, experience 
shows that this iteration only has a minor effect on the 
final stellar parameters and in particular the metallicity, 
as demonstrated in Figure [TJ 

4. METALLICITIES OF STARS WITH AND WITHOUT 
DETECTED TRANSITING PLANETS 

We have determined the metallicities of a sample of 
stars in the Kepler field that have no detected transiting 
planets (SNTP). Although the transit method has not 



Figure 3. The distribution of metallicities of dwarf stars with 
no detected transiting planets (blue, log(g) > 4.2, 88 stars) and 
the sample of stars with small planets (red, 157 stars) from 
IBuchhave et al.1 (|2014T ). The vertical lines show the average metal¬ 
licity of the two samples: [m/H] SNTp dwar f = —0.02 d= 0.02 dex 
and [m/H]g T p = —0.02 d= 0.02 dex, where the uncertainty is the 
standard error of the samples. 

detected any planets, many of the stars in the sample are 
most likely planetary hosts, since the occu r rence rates of 
small planets is high (Howar d et al.l I2012t iFressin et al.l 
120131: iPetigura et al.l I2013[ ). Unfortunately, there is cur¬ 
rently no feasible way to construct a “control sample” of 
stars free of small planetary companions, since our cur¬ 
rent methods are either not sensitive enough or cannot 
exclude the presence of planets. Nevertheless, if some of 
the stars in the SNTP sample do not host planets and if 
such stars are on average more metal poor then planet 
hosting stars, we would expect the average metallicity of 
the SNTP sample to be lower than that of a sample of 
stars hosting planets. 

We compare the metallicities of the SNTP stars to 
samples of stars hosting t ransiting planets (STP) smaller 
than R p < 1.7 i?© from IBuchhave et alJ (1201411 . Short 
period planets receive a large amount of a incident flux 
from their host star and could as a result undergo signifi¬ 
cant evaporation of their atmospheres, which is especially 
true for smaller planets (IQwen fe Wull20l3 l. Such plan¬ 
ets could have accreted a gaseous atmosphere when they 
formed and would later have this atmosphere stripped off 
by photo evaporation. This will decrease the planet radii 
and could therefore potentially distort the metallicity- 
rad ius correlation w e are at tempting to quantify. Follow¬ 
ing IBuchhave et alJ (I2014H , we therefore remove highly 
irradiated planets ( F > 5 x 10 5 J s -1 nW 2 ) from the 
sample. 

Figure [2] shows the distribution of metallicities of the 
two samples (stars with and without detected transiting 
planets). The SNTP sample contains 518 stars and the 
STP sample contains 157 stars. The average metallic¬ 
ity of the two samples is [m/H] gNTP = —0.02 ± 0.01 dex 
and [m/H] STP = —0.02 ± 0.02 dex, where the uncer¬ 
tainty is expressed as the standard error of the mean. 
As mentioned, the stars in the SNTP sample were se¬ 
lected because they have asteroseismic detections. Since 
asteroseismic oscillations are easier to detect in evolved 
stars, the SNTP sample is dominated by sub-giants, 
while the STP sample primarily contains dwarf stars. 
To check whether the dominating sub-giants affect the 
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average metallicity, we exclude stars with weak gravities 
(log(g) < 4.2) from the SNTP sample. Figure [3] shows 
the distribution of the metallicities of the SNTPdwarf 
stars (88 stars) compared with the sample of STP stars. 
The average metallicity of the SNTP sample containing 
only dwarfs is essentially unaffected: [m/H] SNTp dwarf = 
—0.02 I 0.02 dex. 

We performed a two-sample Kolmogorov-Smirnov (K- 
S) test to determine whether the metallicities of the two 
distributions of host stars are not drawn randomly from 
the same parent population. We find a p-value of 0.50 
(0.67 a) when using all the stars in the SNTP sample 
and a p-value of 0.68 (0.41 a) when only including the 
dwarfs. This suggests that the K-S test fails to reject the 
null hypothesis that the two samples are drawn from the 
same parent population. 

We also compare the metallicities of stars hosting 
larger “gas-d warf” planets 1. 7 Rm < R P < 4.0 i?® 
(STPgd) from I Buchhave et al.! ( 2014 ) to the stars with 
no detected transiting planet sample (SNTPdwarf)- The 
average metallicity of the STP g d sample is [m/H] STPgd = 
0.05 ±0.01 dex. A K-S test of the STP gd and SNTP d warf 
samples yields a p-value of 0.0026 indicating that the 
two samples are not drawn from the same parent popu¬ 
lation at a 3.01 a confidence level, in contrast to the STP 
sample of smaller planets with an average metallicity of 
[m/H] STp = —T02 ± 0.02 dex ( in agreement with the 
findings of i Buchhave et al] (1201411 b 

The sample of stars with transiting planets is con¬ 
structed by adding the host star of each individual planet 
to the sample. For multi-planet systems, this can re¬ 
sult in a star being added multiple times if it hosts sev¬ 
eral small planets. If we instead construct a sample of 
unique host stars, with no regard to the number of small 
planets they host, the number of STP stars reduces to 
119 stars. The average metallicity of the STP un j q sample 
is [m/H] STP uniq = ±0.01 ± 0.02 and the K-S p-value is 
0.40 (0.84 a) when comparing to the SNTPdwarf sample. 
Like the other results, the average metallicities of the two 
samples are consistent within the uncertainties and the 
p-value suggests that the K-S test fails to reject the null 
hypothesis that the two sample are drawn from the same 
parent population of stars. 

5. SUMMARY AND DISCUSSION 

We have compared the metallicities of a sample of stars 
with no detected transiting planets (SNTP) to the metal¬ 
licities of a sample of stars hosting sma ll planets (STP 
with R p < 1.7 f?®) from IBuchhave et~ ( 2014* 1. If some 
of the stars in the SNTP sample do not host any planets 
and if these stars are, in fact, more metal poor as a uni¬ 
versal metallicity correlation would suggest, we would 
expect the average metallicity of the sample of SNTP 
stars to be more metal poor than a sample of stars host¬ 
ing small transiting exoplanets. We find the average 
metallicity of the two samples to be [nu/H] g NTP dwarf = 
—0.02 ± 0.02 dex and [m/H] gTP = —0.02 ± 0.02 dex. 

These results seem a t odds with the conclusions from 
I Wang fc Fischeil (1 2015 ). who suggest the existence of a 
universal planet-metallicity correlation extend ing all the 
way d own to the terrestrial planets. I Wang fe Fischer! 
(|2015f ) find a p-value of 0.005 when comparing the metal¬ 
licities of the sample of stars hosting small planets from 


IBuchhave et al.1 (|2014D to the large sample of Kepler In¬ 
put Catalog stars (K1C: [Brown et ahll2011fl with no tran¬ 
siting planet detections, suggesting that the stars in the 
two samples are not drawn from the same parent popu¬ 
lation at a 2.8 a confidence level. The p-value is vastly 
different from our results (p-value of 0.68), which sug¬ 
gests that the samples are, in fact, drawn from the same 
parent population. 

The discrepancy co uld possibly be due to the fact that 
I Wang k. Fischeil (1 201 5l ) compa re metallicities deriv ed us¬ 
ing two different techniques. iBrown et all (120 1 11 1 state 
that the KIC metallicities are 0.17 dex more metal poor 
than spectroscopica ll y der ived metallicities using SME 
(jValenti fe Piskunovl Il996i l, which is attributed to the 
Basian prior used to derive the KIC metallicities be- 
ing rather narrowly pea ked around -0.10 dex. Although 
I Wang fc Fischeil ( 2015 1 did attempt to correct for this 
offset, it can be difficult to compare metallicities de¬ 
rived using different methods, especially when looking 
for subtle differences in metallicities. We also note that 
I Wang fe Fischeil (120151) did not remove the highly irra¬ 
diated planets that may not obey the radius-metallicity 
correlation we are studying (see Section [f} . 

We find the average metallicities of the SNTP stars 
and the STP stars are consistent within the uncertain¬ 
ties. Both samples comprise a large number of stars 
that were analyzed in a homogeneous manner using the 
same tools, thus avoiding a potential bias between the 
samples. A two-sided K-S test yields a p-value of 0.68 
(0.41 <t), suggesting a failure to reject the null hypothe¬ 
sis that the two samples are drawn from the same parent 
population. Other previous publications have reached 
simil a r conclusions dSousa et all 12011 ; IBuchhave et al.l 
l2012t iNeves et all 120131: lEverett et al.l 120 1 3l ( . although 
the publications based on planets detected via radial 
velocities contain only few planets with masses below 
M p <10 M®. We conclude that there is no evidence 
for an enhanced metallicity of Kepler stars hosting small 
planets ( R p < 1.7 i?®) when compared to dwarf stars 
in the Kepler field with asteroseismic detections but no 
detected transiting planets. 
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